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We study K~/K+ ratios as a function of centrality (par- 
ticipant nucleon number), transverse mass (rri(), and rapid- 
ity, in heavy-ion collisions at beam energies between 1A and 
2 A GeV. We use the relativistic transport model that in- 
cludes expicitly the strangeness degrees of freedom and con- 
sider two scenarios for kaon properties in dense matter, one 
with and one without medium modifications of their proper- 
ties. In both scenarios, The K~ /K + ratio does not change 
very much with the centrality, while the K/n and K /ir ratios 
increase with increasing centrality. Significant differences are 
predicted, both in magnitudes and shapes, for the mt spec- 
tra and rapidity distributions of K~ /K + ratio. Experimental 
measurement of these ratios, currently under investigation by 
the FOPI, KaoS, E866, and E895 collaborations, will be useful 
in revealing the kaon in-medium properties. 

25.75.Dw, 26.60. +c, 24.10.Lx 



I. INTRODUCTION 

Whether and how hadronic properties, such as their 
masses, widths, and dispersion relations, are modified in 
hot and dense medium is a topic of great current interest. 
Of particular importance are the medium modifications 
of kaon properties, as they are related to both sponta- 
neous and explicit chiral symmetry breaking, and they 
are useful inputs for the study of kaon condensation and 
neutron star properties ll^. Since the pioneering work 
of Kaplan and Nelson | on the possibility of kaon con- 
densation in nuclear matter, a huge amount of theoretical 
effort has been devoted to the study of kaon properties 
in dense matter, using such diversified approaches as chi- 
ral perturbation theory P-|l2|, the Nambu— Jona-Lasinio 
model and SU(3) Walecka-type mean-field model 

@,^5| ■ Although quantitatively results from these dif- 
ferent models are not identical, qualitatively, a consistent 
picture has emerged; namely, in nuclear matter the K + 
feels a weak repulsive potential, whereas the K~ feels a 
strong attractive potential. 

Experimentally, in-medium properties of kaon and an- 
tikaon can be obtained from the analysis of kaon-nucleus 
scattering |l6|,[l7j and kaonic atom data |18| ]. The infor- 
mation so obtained is, unfortunately, restricted to low 
densities. For the study of kaon condensation, densities 
much higher than that accessible by kaonic atoms are 
involved. This can only be obtained by analysing heavy- 
ion collision data on various observables involving kaon 
and antikaon. Measurements of kaon spectra and flow 
have been systematically carried out in heavy-ion colli- 



sions at SIS (1-2 AGeV), AGS (10 AGeV), and SPS (200 
AGeV) energies [|l9j . The analysis of their yields, spectra, 
and in particular collective flow has indeed provided use- 
ful information about kaon properties in dense nuclear 
matter [pfj[-pc|. So far most of the experimental data 
from the FOPI ||-|3l| and KaoS &H1 collaborations 
at SIS/GSI seem to be consistent with the predictions 
from the chiral perturbation theory. 

To put these conclusions on a firmer footing, additional 
observables and experimental data from independent col- 
laborations will certainly be useful. In this paper we 
study the centrality, transverse mass (m t ), and rapidity 
dependence of the K~ /K + ratio in heavy-ion collisions. 
Since the medium effects act oppositely on kaon and an- 
tikaon, their ratio can reflect more precisely these effects, 
as was pointed out in Refs. (|^J3^]. We will consider 
Ni+Ni, Ru+Ru and Au+Au collisions at beam energies 
between 1 and 2A GeV. These systems are being inves- 
tigated by the FOPI (3f| and KaoS |3|] collaborations 
at SIS/GSI. Furthermore, the Au+Au collisions are also 
being analysed by the E866 [fH§7|] and E895 @,|l| col- 
laborations at the AGS/BNL. 

This paper is arranged as follows. In Section II, we 
briefly review the relativistic transport model, kaon in- 
medium properties and elementary kaon production cross 
sections. The results are presented in Section III. The 
paper ends with a short summary in Section IV. 



II. THE RELATIVISTIC TRANSPORT MODEL 
AND KAON PRODUCTION 

Heavy-ion collisions involve very complicated nonequi- 
librium dynamics. One needs to use transport models in 
order to extract from experimental data the information 
about in-medium properties of hadrons. In this work we 
will use the relativistic transport model similar to that 
developed in Ref. |^]. Instead of the usual linear and 
non-linear a-u; models, we base our model on the effective 
chiral Lagrangian recently developed by Furstahl, Tang, 
and Serot [Q, which is derived using dimensional anal- 
ysis, naturalness arguments, and provides a very good 
description of nuclear matter and finite nuclei. In the 
mean-field approximation, the energy density for the gen- 
eral case of asymmetric nuclear matter is given by 
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The nucleon effective mass m* N is related to its scalar 
field $ by m* N = thin — W and i? are the isospin-even 
and isospin-odd vector potentials, respectively. The last 
three terms give the self-interactions of the scalar field, 
the vector field, and the coupling between them. The 
meaning and values of various parameters in Eq. (1) can 
be found in @. 

From the energy density of Eq. (1), we can also de- 
rive a relativistic transport model for heavy- ion collisions. 
At SIS energies, the colliding system consists mainly of 
nucleons, delta resonances, and pions. While medium 
effects on pions are neglected, nucleons and delta reso- 
nances propagate in a common mean-field potential ac- 
cording to the Hamilton equation of motion, 
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dp 

~dt 



= -V X {E* + W), 
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where E* = VP* 2 +m* 2 . These particles also undergo 
stochastic two-body collisions, including both elastic and 
inelastic scattering. 

In heavy-ion collisions at incident energies considered 
in this work, kaons can be produced from pion-baryon 
and baryon-baryon collisions. For the former we use cross 
sections obtained in the resonance model by Tsushima 
et al. Q]. For the latter the cross sections obtained in 
the one-boson-exchange model of Ref. |l5],[l6| are used. 
Both models describe the available experimental data 
very well. For antikaon production from pion-baryon col- 
lisions we use the parameterization proposed by Sibirtsev 
et al. Q. For baryon-baryon collisions, we use a some- 
what different parameterization, which describes the ex- 
perimental data better, than Ref. |47f| . In addition, the 
antikaon can also be produced from strangeness-exchange 
processes such as irY — > KN where Y is either a A or £ 
hyperon. The cross sections for these processes are ob- 
tained from the reverse ones, KN — > ttY , by the detailed- 
balance relation. All the parameterizations for the ele- 
mentary cross sections and comparisons with experimen- 
tal data can be found in our recent paper |2q| . 

Particles produced in elementary hadron-hadron inter- 
actions in heavy-ion collisions cannot escape the environ- 
ment freely. Instead, they are subjected to strong final- 
state interactions. For the kaon, because of strangeness 
conservation, its scattering with nucleons at low energies 
is dominated by elastic and pion production processes, 
which do not affect its final yield but changes its mo- 
mentum spectra. The final-state interaction for the an- 
tikaon is much stronger. Antikaons can be destroyed in 
the strangeness-exchange processes. They also undergo 



elastic scattering. Both the elastic and absorption cross 
sections increase rapidly with decreasing antikaon mo- 
menta. This will have strong effects on the final K~ 
momentum spectra in heavy-ion collisions. 

we will consider two scenarios for kaon properties in 
nuclear medium, one with and one without medium mod- 
ification. From the chiral Lagrangian the kaon and an- 
tikaon in-medium energies can be written as 28 



ojk = [m K + k 2 - a K ps + (bKp) 2 ] + b K p (3) 

w iC = [ m K + k 2 - a R p s + (bxp) 2 ] 1/2 - b K p (4) 

where bx — 3/(8/ 2 ) w 0.333 GeVfm 3 , ajf and are 
two parameters that determine the strength of the at- 
tractive scalar potential for kaon and antikaon, respec- 
tively. If one considers only the Kaplan-Nelson term, 
then clk — »if = ^kn/ '/ 2 - In the same order, there is 
also the range term which acts differently on kaon and 
antikaon, and leads to different scalar attractions. Since 
the exact value of and the size of the higher-order 

corrections are still under intensive debate, we take the 
point of view that a K ^ can De treated as free parameters 
and try to constrain them from the experimental observ- 
ables in heavy-ion collisions. In Ref. (2^j2^] we found 
that a K fa 0.22 GeV 2 fm 3 and w 0.45 GeV 2 fm 3 pro- 
vide a good description of kaon and antikaon spectra in 
Ni+Ni collisions at 1A and 1.8A GeV. These values will 
be used in this work as well. 
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FIG. 1. Proton transverse mass spectra in central Ni+Ni 
collisions at 1.93 AGeV in three rapidity bins. The open 
circles are experimental data from the FOPI collaboration 
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In nuclear medium, kaon and antikaon masses are mod- 
ified, so are their production thresholds. We will then use 
a/sJJ, which are calculated with effective masses, in eval- 
uating the in-medium production cross sections. This 
amounts to the change of threshold, or approximately, 
to the change of available phase space. In addition to 
the change in the production cross sections, the medium 
effects on kaon and antikaon also affect their momentum 
spectra, when they propagate in the mean-field poten- 
tials. The Hamilton equations of motion for kaon and 
antikaon are very similar to those for nucleons pl[|, 

dr k dk _ ,_. 

at uj K j i + bkPN at 

where the minus sign corresponds to kaon, and the plus 
sign to antikaon. It is clearly that the K + momentum 
increases and that of K~ decreases when they propagate 
in their respective mean field potentials. This affects sig- 
nificantly their momentum spectra, and especially the 
momentum spectra of their ratio. 
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FIG. 2. Same as Fig. 1, for 7T~. 



III. RESULTS AND DISCUSSIONS 

In this section we present our results for the centrality, 
transverse mass, and rapidity dependences of K~ / K + 
ratio. Before that, we first compare the proton and pion 
mt spectra from our calculation with the available exper- 
imental data, taking central Ni+Ni collisions at 1.93 A 
GeV as an example. In Fig. 1 we compare the proton 
transverse mass spectra in three rapidity bins with the 



experimental data from the FOPI collaboration p^Jf9| |. 
Similar comparison for the 7r~ transverse mass spectra 
is shown in Fig. 2. Our results are seen to be in good 
agreement with the data Similarly, in Ref. |2|] 

the proton and 7r~ rapidity distributions obtained in our 
transport model were shown to be in good agreement 
with the FOPI data. 



A. centrality dependence 

In this subsection, we discuss the centrality depen- 
dence of the K~ /K + ratio. We use the participant nu- 
cleon number A part as the measurement of the centrality. 
In Fig. 3 we show the A part dependence of the pion mul- 
tiplicity in Ni+Ni collisions at 1.93A GeV. The solid line 
gives 1 /3 of the total pion number obtained in our calcu- 
lation, while the circles and squares are the FOPI data 
|5(if for 7r~ and ir + , respectively. It is seen that our re- 
sults are in good agreement with the data, and that the 
pion multiplicity increases almost linearly with the par- 
ticipant nucleon number. In other words, the N^/Apart 
ratio is almost independent of the centrality. 
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FIG. 3. Pion multiplcity as a function of the participant 
nucleon number in Ni+Ni collisions at 1.93 AGeV. 



On the other hand, it is well-known that the K + (usu- 
ally shown in terms of K + /ir + ratio) yield increases more 
than linearly with the A part for beam energies ranging 
from 1A GeV H-HH to 10A GeV US ° U1 ' results 
for K/ir and K/n ratios in Au+Au collisions are shown 
in Fig. 4. It is seen that either with or without kaon 
medium effects, these ratios increase more than linearly 
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with the participant nucleon number. This is due to the 
increasing importance of the secondary processes involv- 
ing baryon-resonances, pions, and hyperons in the case 
of antikaon production, when going from peripheral to 
central collisions. The E866 data for Au+Au collisions 
at 11. 6 A GeV/c indicate that K~ /ir + ratio also increases 
with increasing centrality ^jj . 

4 | , 1 , 1 , 1 , 1 

Au+Au, 2A GeV 

no kaon medium effects 

3 — with kaon medium effects 

I 

2 - - 




0.2 - 



0.1 - 



100 200 300 400 

"part 



5 



4 







1 


1 1 1 1 1 1 

Ni + Ni, 1.8A GeV 


1 1 1 


- , 






c 


) 


- 




C 


> 






( 


> 






( 


3 








- 


no kaon 


medium effects 




with kaon medium effects 




O KaoS 






1 




, 1 




1 



I I I I I I I I I I I I I 

20 40 60 80 100 120 



part 

FIG. 5. K~ /K + ratio as a function of the participant nu- 
cleon number in Ni+Ni collisions at 1.8A GeV. The open cir- 
cles are the experimental data from the KaoS collaboration 
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FIG. 4. K/tt and K/tt ratios as a function of the partici- 
pant nucleon number in Au+Au collisions at 2A GeV. 



The centrality dependence of the K~/K + ratios is 
shown in Figs. 5, 6, and 7, for Ni+Ni at 1.8A GeV, 
Ru+Ru at 1.69A GeV, and Au+Au at 2A GeV, respec- 
tively. For Ni+Ni collisions, we also shown in Fig. 5 the 
experimental data from the KaoS collaboration jj5| by 
open circles. For all the systems considered, it is seen 
that either with or without medium effects, the K~ /K + 
ratio does not depend very much on the centrality, since 
in both scenarios, the K/tt and K/tt ratios increase at 
about the same rate with increasing centrality (see Fig. 
4) . The experimental data of the KaoS collaboration for 
Ni+Ni collisions (Fig. 5) show little increase of the ratio 
towards central collisions within their statistical uncer- 
tainties. Preliminary data from the FOPI collaboration 
for Ru+Ru at 1.69 A GeV Q, and data from E866 col- 
laboration for Au+Au collisions at 10. 6A GeV |36| both 
show weak centrality dependence for the K~ / K + ratio. 
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FIG. 6. Same as Fig. 5, for Ru+Ru collisions at 1.69A 
GeV. 



Without kaon medium effects, the K / K + ratio 
ranges from 0.006 in Ru+Ru at 1.69A GeV, to 0.01 
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in Ni+Ni collisions at 1.8A GeV, to 0.0075 in Au+Au 
collisions at 2A GeV. The increase of this ratio from 
Ru+Ru to Ni+Ni comes from the increase in beam en- 
ergy, since at these energies the antikaon excitation func- 
tion is steeper than that of kaon because of a higher 
threshold. On the other, the decrease of this ratio from 
Ni+Ni to Au+Au is due to the fact that in the large 
Au+Au system the antikaon absorption effects become 
more significant. Qualitatively, the absorption proba- 
bility is proportional to the product of the absorption 
cross section and the average size of the system. In this 
case the experimental data from the KaoS collaboration 
for Ni+Ni collisions are significantly underestimated (see 
Fig. 5). 
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FIG. 7. Same as Fig. 5, for Au+Au collisions at 2A GeV. 



When kaon medium effects are included, the K~ / K + 
ratio increase by about a factor of 4, for all the systems 
considered here. This increase results from a factor of 
3 increase in antikaon yield and about 35% reduction in 
kaon yield. In this case, the K~ /K + ratio from the KaoS 
collaboration can be nicely explained. Furthermore, pre- 
liminary data from the FOPI collaboration indicate a 
K~ /K + ratio of about 0.02-0.03 in Ru+Ru collisions at 
1.69 A GeV. Our predictions that include kaon medium 
effects are seen to be in better agreement with these pre- 
liminary data than those without the kaon medium ef- 
fects. 

Naively, it is expected that when kaon medium effects 
are included, the K~ /K + ratio should increase with the 
increasing centrality, since in central collisions reduction 
of the antikaon mass and the increase of kaon mass are 
the most significant. However, since the second processes 



Yn — » KN play an important role in antikaon produc- 
tion, and hyperon yield, which are produced in associa- 
tion with kaons, is reduced most significantly in central 
collisions, the increase in the reduction of antikaon mass 
towards central collisions is largely compensated. The 
increase of the K~ /K + ratio with increasing centrality 
when kaon medium effects are included is thus marginal 
(on the order of 10-20%). 

It is also of interest to show the beam energy depen- 
dence of the K" / K + ratio in central Au+Au collisions. 
This is done in Fig. 8 for impact parameter b= 1 fm. In 
both scenarios with and without kaon medium effects, the 
ratio is seen to increase as beam energy increases. This 
is understandable, as the antikaon production threshould 
is higher than that of kaon, so at these energies, the 
antikaon production cross sections increases faster than 
that of kaon. 

Finally, we show in Fig. 9 the time evolution of cen- 
tral density, kaon yield, and antikaon yield in Au+Au 
collisions at 1.5A GeV and b=3 fm. It is seen that for a 
considerable duration of time the system is compressed 
to a density of around 3po, and during this period of 
time, most of the kaons and antikaons are produced. As 
was shown in Ref. [^8| , K~ condensation is predicted to 
occur around three times normal nuclear matter based 
on current theoretical and empirical information. The 
experimental data on kaon and antikaon yields, spectra, 
ratios, and flow from Au+Au collisions at 1.5A Gev will 
be very useful in pinning down the critical density for 
kaon condensation. 
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FIG. 8. Beam energy dependence of K / K + ratio in cen- 
tral Au+Au collisions. 
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FIG. 9. Time evolution of central density as well as kaon 
and antikaon yields in Au+Au collisions at 1.5A GeV and 3 
fm. 



momenta, leading to an enhanced K~ /K + ratio at small 
transverse masses. From Fig. 10, we also see that the 
effects of propagation in mean-field potentials are more 
pronounced in central collisions. 
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B. Transverse mass spectra 

The effects of kaon and antikaon mean-field potentials 
can be more clearly seen by looking at their ratio as a 
function of their tranverse mass (or kinetic energy). The 
results are shown in Figs. 10, 11, 12, and 13, for Ni+Ni 
at 1.8A GeV, Ru+Ru at 1.69A GeV, Au+Au at 1.5A 
GeV, and Au+Au at 2A GeV, respectively. The results 
for Ni+Ni are presented in terms of kinetic energy for 
minimum-biased collisions, in accordance with the ex- 
perimental data from the KaoS collaboration |}5j , shown 
in Fig. 10 with open circles. 

When kaon medium effects are neglected, the K~ / K + 
ratio is seen to increase slightly with increasing transverse 
momentum. Since the antikaon absorption cross sec- 
tion by nucleons becomes large at low momentum, low- 
momentum antikaons are more strongly absorbed than 
high- momentum ones. This makes the K~ / K + ratio de- 
crease at small transverse mass. When medium effects 
are included, we find that the shape of the K~ / K + ratio 
is completely different from that without kaon medium 
effects. The ratio is now seen to increase dramatically to- 
wards small transverse mass. For example, for Ru+RU 
collisions, it increases from about 0.02 at m t — vhk = 0.3 
GeV to about 0.08 when m t — vtlk approaches zero. The 
difference in the shape of the m t spectra comes from the 
propagation of kaons and antikaons in their mean-field 
potential. Kaons are 'pushed' to high momenta by the 
repulsive potential, while antikaons are 'pulled' to low 



FIG. 10. Kinetic energy spectra of the K~ / K + ratio in 
Ni+Ni collisions at 1.8A GeV. The experimental data are from 
the KaoS collaboration [B5l 
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FIG. 11. Transverse mass spectra of K /K + ratio in cen- 
tral Ru+Ru collisions at 1.69A GeV. 
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FIG. 12. Same as Fig. 11, for Au+Au collisions at 1.5A 
GeV. 



Experimental data from the KaoS collaboration pro- 
vided some indication for this dramatic increase of the 
K~/K + ratio at small momenta ]35|j . Our results are 
compared with these data in Fig. 10. It is seen that our 
results including the kaon medium effects are in much 
better agreement with the data. Of course the statistical 
uncertainty of the data is still quite large (especially for 
the first datum point at 0.12 GeV). This will be improved 
in their recent analysis of Ni+Ni collisions at 1.93 A GeV 
p9| . It will also be very useful if the ratio at Identic ener- 
gies less than 0.1 GeV can be measured experimentally. 

It should be emphasized that in this work, as in our 
previous study of kaon production in heavy-ion collisions, 
the explicit momentum dependence of the kaon scalar 
and vector potential is not considered. This should be 
a reasonable approximation for heavy-ion collisions at 1- 
2A GeV, as kaons (antikaon) produced in these reactions 
usually have small momenta. This approximation will 
break down at higher beam energies about 10A GeV. A 
recent theoretical calculation |5^] indicated that the at- 
tractive antikaon potential becomes weaker as its momen- 
tum (relative to the medium) increases. Experimental 
data on ^"-nucleus scattering at high incident momenta 
also indicate a repulsive (rather than an attractive) an- 
tikaon optical-model potential. This momentum depen- 
dence of the kaon (antikaon) potential might explain the 
fact that so far the experimental data for heavy-ion colli- 
sions at 10A GeV do not show substantial enhancement 
of the K~ / K + ratio approaching small transverse masses 
(momenta) . 
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FIG. 13. Same as Fig. 11, for Au+Au collisions at 2A GeV. 



C. Rapidity distribution 

The rapidity distributions of the K~ /K + ratio should 
provide quite similar information on kaon medium effects 
as its transverse mass spectra. The results for Ni+Ni at 
1.93A Gev, Ru+Ru at 1.69A GeV, Au+Au at 1.5A GeV, 
and Au+Au collisions are shown in Figs. 14, 15, 16, and 
17, respectively. For Ni+Ni collisions we show also the 
preliminary data from the FOPI collaboration |]52] , p4]| by 
circles. 

Indeed, the magnitudes and the shapes of these rapid- 
ity distributions with and without kaon medium effects 
are very different. Without kaon medium effects, the ra- 
pidity distribution is seen to decrease slightly from tar- 
get/projectile rapidities to mid-rapidity. This is again 
due to the large absorption cross section for slow-moving 
antikaons. When kaon medium effects are included, there 
appears not only an overall increase of the ratio, because 
of the increased production cross section of antikaons, 
but also a significant change of the shape of the rapid- 
ity distribution. The K~ / K + ratio is seen to increase 
steadily from target /projectile rapidities to mid-rapidity. 
Our results with kaon medium effects are in much bet- 
ter agreement with the preliminary data from the FOPI 
collaboration (Fig. 14). 
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FIG. 14. Rapidity distribution of K~/K+ ratio in central 
Ni+Ni collisions at 1.93A GeV. The solid circles are prelimi- 
nary experimental data from the FOPI collaboration 
while the open circles are obtained by reflecting the data with 
respect to the mid-rapidity. 
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FIG. 16. Same as Fig. 14, for Au+Au collisions at 1.5A 
GeV. 



For Au+Au collisions at 2A GeV and 6 = 1 fm, we also 
show in the figure the K~ /K + ratio before any final-state 
interactions on kaons and antikaons are included, which 
means that even the antikaon absorption is turned off. 
We see that in this case, the ratio also decreases from 
mid-rapidity to projectile and target rapidities. The de- 
crease of this ratio towards large rapidities (momenta) 
reflects the fact that antikaons are produced at a higher 
threshold than kaons, so that their momentum distribu- 
tions are restricted more severely by the available en- 
ergies. However, we need to emphasize that antikaon 
absorption and the fact that its absorption cross section 
increases at low momenta are well-known experimental 
facts and must be included in transport model calcula- 
tions. 

As mentioned in the Introduction, so far most of the 
experimental data from the FOPI and KaoS collabora- 
tions at SIS/GSI are consistent with the chiral pertur- 
bation theory predictions for kaon in-medium properties. 
For an ultimate confirmation of these medium effects, 
it should be be very useful if independent experimental 
data for Au+Au collisions at 2A GeV from the E866 and 
E895 collaborations could become available. 



FIG. 15. Same as Fig. 14, for Ru+Ru collisions at 1.69A 
GeV. 
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FIG. 17. Same as Fig. 14, for Au+Au collisions at 2A GeV. 



IV. SUMMARY 

In summary, we studied K~ /K + ratios as a function of 
centrality (participant nucleon number), transverse mass 
(m t ), and rapidity, in Ni+Ni, Ru+Ru, and Au+Au colli- 
sions at beam energies between 1 and 2A GeV. We used 
the relativistic transport model that includes explicitly 
the strangeness degrees of freedom and considered two 
scenarios for kaon properties in dense matter, one with 
and one without medium modifications of their prop- 
erties. In both scenarios, The K~/K + ratio does not 
change very much with the centrality, while the K/tt and 
K/tt ratios increase with increasing centrality. Signifi- 
cant differences were predicted, both in magnitudes and 
shapes, for the m t spectra and rapidity distributions of 
K~/K + ratio. We found that the experimental data 
from the KaoS collaboration for the kinetic energy spec- 
tra of the K~ /K + ratio and those from the FOPI col- 
laboration for its rapidity distribution support the sug- 
gestion of kaon medium effects. Wc emphasize that the 
independent data from the E866 and E895 collaborations 
for Au+Au collisions at 2A GeV will be useful in confirm- 
ing or confronting these findings. 

We are grateful to N. Herrmann and P. Senger for send- 
ing us the data files, and to N. Herrmann, C. Ogilvie, and 
P. Senger for useful communications. This work is sup- 
ported in part by the Department of Energy under Grant 
No. DE-FG02-88ER40388. 
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